:s«SKHr?sar — 


<T<rJU*^ 


1  ^^fcajfcod* 


mm 


: 


';  *ei  x  no 


k-V'}^  'r, &*■■  ■  ■  ■'■  •• 


Ref 


<  ”*Ut 


J  fcJ- 

i 


Reft1  -<»  TWftt  ^tftFUCijJ  W|  -5«  00;z  '•«  •  *  -  • 

4^SS».w- 

• ‘::  ;  >"?4 


)  Abstract  .onfyr^ot  puW.liHeti 
B6pX*chap(*r 
Cbfforsnc&^Jbceed  ingt 
(hbtfefented). 


STftN  imOlA/rert*  OS-  9039 
Rovte.Stoet  No.7tegy 


Job OrtietfNo  J3r5p2^E8ji5  .  _ 


This  is  aHr^i  Security  Review 
Any Ranges  made  iivthe  document 
after  approved  by  Cede  1 2*.6 
nullify  the  Security  Review 


,  \  ’.  \titimeof  Canfcrencrf 

'■•  ;• .  '■■  '■  ;  *;;f  '  rebate,  P&tre  EST Ctasstno&tlon  of  Conference; 


and/or  tar'^iaetjon:^  Oca^ra>1tf<rtinrt^t  Ufleiaritfflod 


:  W&tiM  C/^iriifea  fofi  ;  of-  PutifkaMn) . 


rw»rtre,W-P06«8»»''J. 

After pT^*Bffrtfertfor  qr  puWlc^fiofl,  pe^tfbrrtpOtHfc^  will  be  -entered  -In  the  publicattohis  database,  In  accoitJfence 

!tte;lh€^  _ )  (li  ho* _ *)  dasSlfled.t^  refB^ni^tb). 

TW^p'^lp^T^^^^TTOtvtfeteta^t?y^isdfesii^pCb§^dft^^^al^-or*l»ugg^dlibTO^f  outsWe-ln'dividuifls-  or  conbem^wWdi-baye  beeh 
.cdffiifiwi&aSiyto  paber^does  V:».  Yfdoerndtv  XVcbntah'onv^infibflfe^^ite 

\1  v-  IfhA-faver  ^&.v -Vbde n facQrbdfatbd tfr  amvdffi^ Report^ 


-  *  ■  jjl. 


N#«dbv 


^u&llcly  totfrcw  darter  ttvsfri  button 


LiLrJl^V 


jrpSoM^tfWHj4p*t  n 

*  B-*,»>-  t"-  V  -  \!  T  - 

T  p  iWtff  Mil  ICFPW'TnJy  ** 


ffiSwfld 


ofm*  Stanton,  lb* 


fcbfajw*  mbn»f  orti^PWKfr  (b«  -. _ ) 

ffa*  nw>  ^  1  ttoOn'd*»H!«L 


t  M-lM 


\;.Pqrt  or*&SVit**T\%  mi&Mfcti.  .  r./*\  ry 

•  Acdpf^uJttW  toifcvbftc*.  */*£?*  0 


mu 


z- 


-  .-» ■  "O 


ftW  •  -  ; ' ;  . .  ;  %:  ■  .;•  'v 


ELSEVIER 


Available  online  at  www.sciencedirect.com 

ScienceDirect 


Ocean 

Modelling 


Ocean  Modelling  22  (2008)  48  -65  _ _ 

www.clscvicr.com/locatc/occmod 


Relative  dispersion  from  a  high-resolution  coastal  model 

of  the  Adriatic  Sea 

Angelique  C.  Hazaa  *,  Andrew  C.  Poje,  Tamay  M.  Ozgokmen3,  Paul  Martin0 

a  RSMASIMPO.  University  of  Miami,  Miami.  FL,  USA 
b  Department  of  Mathematics ,  College  of  Staten  Island,  CUN  Y.  Staten  Island,  USA 
c  Naval  Research  Laboratory.  Stennis  Space  Center,  MS.  USA 

Received  22  August  2007;  received  in  revised  form  10  January  2008;  accepted  1 1  January  2008 
Available  online  6  February  2008 


Abstract 

Synthetic  drifter  trajectories  computed  from  velocity  data  produced  by  a  high-resolution  NCOM  model  are  used  to  investigate  the 
scaling  of  relative  dispersion  and  the  distribution  of  finite-scale  Lyapunov  exponent  (F'SLE)  fields  in  the  Adriatic  Sea.  The  effects  of  vary¬ 
ing  degrees  of  spatial  and  temporal  filtering  of  the  input  Eulerian  velocity  fields  on  the  Lagrangian  statistics  are  investigated  in  order  to 
assess  the  sensitivity  of  such  statistics  to  model  error.  It  is  shown  that  the  relative  dispersion  in  the  model  Adriatic  circulation  is  generally 
super-diffusive,  sealing  nearly  ballistieally  in  elose  agreement  with  Lagrangian  observations  from  a  limited  set  of  drifters.  The  large-scale 
dispersion  is  dominated  by  persistent  separation  regions  and  the  controlling  influence  of  the  Western  Adriatic  Current  (WAC).  Temporal 
filtering  with  averaging  windows  up  to  monthly  time  scales  only  affects  the  relative  dispersion  at  scales  smaller  than  20  km  without  alter¬ 
ing  the  overall  scaling  regime.  In  contrast,  spatial  smoothing  at  scales  as  small  as  5  km  significantly  reduces  relative  dispersion  at  all  scales 
up  to  100  km.  While  basin-scale  dispersion  statistics  are  strongly  dependent  on  spatial  resolution  of  the  model  WAC,  maps  of  FSLE 
fields  over  initial  conditions  indicate  that  the  detailed  geometry  of  the  dispersion  is  determined  to  a  large  extent  by  the  temporal  reso¬ 
lution  of  the  model.  In  addition,  the  degree  of  spatial  heterogeneity  in  the  flow  field  implies  that  the  existence,  or  non-existence,  of  a 
distinct  exponential  regime  in  the  FSLE  at  small  scales  is  extremely  sensitive  to  the  details  of  particle  pair  sampling  strategies. 
Published  by  Elsevier  Ltd. 


1.  Introduction 

The  transport,  spreading,  and  mixing  of  passive  particles 
and  tracers  have  a  wide  range  of  practical  applications  in 
both  classical  and  geophysical  fluid  flows.  These  include 
the  spreading  of  industrial  plumes  in  the  air,  of  pollutants, 
discharges,  and  fish  larvae  in  the  ocean,  or  global  adveetion 
of  temperature  and  tracers  in  the  ocean  thermohaline  circu¬ 
lation.  The  Lagrangian  framework 

^  =  v(0=u(x,0  (1) 

is  the  natural  setting  for  such  problems,  where  x  is  the  po¬ 
sition,  v  is  the  velocity  of  the  Lagrangian  particle,  and  u  is 
the  velocity  field  of  the  underlying  Eulerian  flow. 
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Many  of  the  above-mentioned  practical  applications  can 
be  thought  of  in  terms  of  the  adveetion  of  a  cluster  of 
Lagrangian  particles  with  the  flow  field.  One  of  the  metrics 
that  can  be  used  to  describe  this  problem  is  the  absolute  (or 
single-particle)  dispersion 

A 2('-'o)  M)  -  =  (!*'(')  -  *.('o)|2>,  (2) 

1  =  1 

which  quantifies  the  average  distance  covered  in  time  by  N 
Lagrangian  particles  following  their  release  at  t  =  /0. 

Another  Lagrangian  metric  more  closely  tied  to  scalar 
mixing  processes  is  the  relative  (two-particle)  dispersion 

D2(0  =  (|x'1>(/)-x‘2>(/)|2).  (3) 

where  the  average  is  over  all  pairs  of  particles  in  the  cluster. 
Thus,  relative  dispersion  quantifies  the  mean-square  parti¬ 
cle  separation  and  is  the  focus  of  the  present  study. 
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The  interpretation  of  relative  dispersion  defined  in  (3)  is 
often  obseured  in  the  presence  of  intermittent  turbulenee 
and  finite-sized  domains  (BofTetta  et  al.,  2000).  The  root 
eause  of  this  is  that  the  average  in  (3)  is  taken  over  a  wide 
range  of  separation  seales  at  any  one  time  and  is  thus  dom¬ 
inated  by  those  particle  pairs  with  the  largest  separations. 
To  address  this,  a  complementary  measure,  namely  the 
finite-scale  Lyapunov  exponent  (FSLE),  is  often  employed 
(Artale  et  al.,  1997;  Aurell  et  al.,  1997) 


A(<5)  = 


1n(g) 

<*(«)>’ 


(4) 


where  (i(<5))  is  the  averaged  time  (over  the  number  of  par- 
tiele  pairs)  required  to  separate  from  a  distance  of  S  to  xS. 
Like  Richardson’s  original  time-distanee  graphs  (Riehard- 
son,  1926),  the  FSLE  provides  a  measure  of  dispersion  as  a 
function  of  the  spatial  seale  and  serves  to  isolate  the  differ¬ 
ent  dispersion  regimes  corresponding  to  the  different  seales 
of  the  oeeanie  flows.  Spatial  distributions  of  the  FSLE  and 
related  measures  sueh  as  finite-time  Lyapunov  exponents 
have  recently  been  used  to  define  so-ealled  Lagrangian 
coherent  structures  (Shadden  et  al.,  2005)  and  to  map  trans¬ 
port  barriers  and  mixing  regimes  in  a  number  of  oeeanie 
applications  (Abraham  and  Bowen,  2002;  d’Ovidio  et  al., 
2004;  Waugh  et  al.,  2006;  Olaseoaga  et  al.,  2006). 

The  study  of  relative  dispersion  in  fluid  flows  is  a  funda¬ 
mental  area  of  research,  dating  baek  to  the  1920s.  The 
reader  is  referred  to  Bennett  (1987)  and  Sawford  (2001) 
for  a  comprehensive  review  of  the  subjeet  matter.  Here,  a 
short  summary  of  this  material  is  provided  followed  by 
how  the  present  study  is  motivated  and  why  it  is  original 
in  the  eontext  of  oeeanie  flows. 

Traditionally,  relative  dispersion  has  been  related  to  tur¬ 
bulent  easeade  processes  in  the  Eulerian  velocity  field  (Ben¬ 
nett,  1984;  Piterbarg,  2005).  The  original  formula  proposed 
by  Riehardson  (1926),  D2(t)  =  where  C  =  0(  1)  is  a 
constant  and  t  is  the  turbulent  kinetie  energy  dissipation 
rate,  is  found  to  be  consistent  with  Kolmogorov’s  (1941) 
energy-easeade  law  for  three-dimensional  (3D),  homoge¬ 
neous,  isotronie  turbulenee  in  the  inertial  range, 
E(k )  ~  t2hk  (Batehelor,  1952).  Nevertheless,  the  precise 
mechanisms  responsible  for  Richardson’s  (1926)  relation 
and  the  range  of  conditions  under  which  it  is  valid  remains 
the  subjeet  of  some  debate  even  in  idealized  experimental 
and  numerical  settings  (e.g.,  Bourgoin  et  al.,  2006).  The 
direet  applicability  to  oeeanie  flows,  whieh  tend  to  be 
two-dimensional  (2D)  under  the  effects  of  rotation  and 
stratification  (Cushman-Roisin,  1995)  and  whieh  are  often 
dominated  by  persistent  inhomogeneity  and  anisotropy,  is 
open  to  question.  Kraiehnan  (1967)  proposed  that  the  sim¬ 
plest  2D  turbulence  ean  exhibit  a  dual  spectrum,  namely  an 
enstrophy  easeade  E(k)  f'k  3,  where  /?  is  the  enstrophy 
dissipation  at  seales  smaller  than  the  foreing  seale,  k  >  /cf, 
and  an  inverse  energy-easeade  E(k )  e2nk~5/  1  at  seales  lar¬ 
ger  than  the  foreing  seale,  k  <  kr.  In  the  enstrophy  easeade 
regime,  exponential  partiele  separation  ean  be  expeeted, 
Er(t)  —  D~  exp(/o/),  where  is  the  enstrophy  easeade 


time  seale  (Lin,  1972).  Note  that  the  emergence  of  a  dual 
speetrum  in  2D  turbulenee  requires  foreing  at  a  narrow 
range  of  seales  and  the  speetrum  tends  toward  E(k)  ~  k  3 
for  the  more  typieal  oeeanie  ease  of  broad-band  foreing 
(Lesieur,  1997).  At  the  asymptotie  limit  of  very  long  times 
and  partiele  separations,  when  all  partiele  correlations 
diminish,  the  so-ealled  diffusive  regime,  D2(t)  =  2k7,  where 
k  is  the  diffusivity  eoeffieient,  should  be  attained  (Taylor, 
1921). 

Reeent  developments  in  the  understanding  of  Lagrang- 
ian  transport  and  mixing  (Wiggins,  2005)  indicate  that 
there  are  issues  other  than  turbulent  easeades  that  deserve 
consideration  for  understanding  oeeanie  relative  disper¬ 
sion  Of  particular  importance  is  the  eoneept  of  ehaotie 
adveetion,  in  whieh  (laminar)  Eulerian  flow  fields  with  rel¬ 
atively  simple  spatial  structure  ean  lead  to  eomplex 
Lagrangian  behavior  due  to  the  effeet  of  time  dependence 
alone  (Aref,  1984).  Chaotie  flows,  by  their  very  definition, 
lead  to  exponential  partiele  separations  for  small  initial 
displacements. 

Finally,  various  empirieal  correlations  have  been  estab¬ 
lished  between  the  Lagrangian  relative  dispersion  and  the 
type  of  coherent  structures  in  the  Eulerian  fields.  For 
instance,  the  ballistic  regime,  D2(t)  ~  t2,  is  observed  in  the 
presence  of  persistent  correlations  and  shear  dispersion 
(Iudieone  et  al.,  2002)  and  the  sub-diffusive  regime, 
D2{t)  ty  with  y  <  1,  is  often  attributed  to  partiele  trapping 
within  eddies  (Cardoso  et  al.,  1996;  Provenzale,  1999). 

While  it  is  possible  to  parse  further  the  spatial  and  tem¬ 
poral  seales  into  other  regimes  of  dispersion  on  a  theoreti¬ 
cal  basis  (Sawford,  2001;  Piterbarg,  2005),  there  appears  to 
be  a  high  degree  of  uncertainty  about  whieh  regimes  are 
actually  attained  in  praetiee,  namely  in  the  laboratory, 
atmospheric  and  oeeanie  observations,  and  numerical  sim¬ 
ulations.  Eaeh  one  of  these  methods  introduces  various  fac¬ 
tors  that  ean  eause  deviations  from  theory.  For  instance, 
the  Reynolds  number,  thus  the  range  of  the  interacting  tur¬ 
bulent  processes,  tends  to  be  limited  in  laboratory  experi¬ 
ments  and  direet  numerical  simulations  in  comparison  to 
that  in  the  real  oeean.  Also,  the  ocean  is  subjeet  to  a  very 
complex  foreing  at  the  surfaee,  the  presence  of  strongly 
inhomogeneous  mean  currents,  and  the  limiting  effeets  of 
the  boundaries.  These  effeets  are  expeeted  to  lead  to  devia¬ 
tions  of  relative  dispersion  from  the  regimes  that  are 
expeeted  on  the  basis  of  theory.  As  sueh,  a  brief  review 
of  the  observational  experience  follows. 

In  laboratory  experiments  of  2D  homogeneous  turbu¬ 
lenee,  Julien  (2003)  reported  a  elear  exponential  separation 
regime.  Bourgoin  et  al.  (2006)  obtained  Batehelor  (1952) 
sealing  (similar  to  ballistic)  in  3D  homogeneous  turbulenee 
in  the  laboratory,  but  eould  not  eonfirm  the  Riehardson 
regime.  Numerical  simulations  of  2D  homogeneous  turbu¬ 
lenee  (Boffetta  and  Celani,  2000;  Braeeo  et  al.,  2004) 
yielded  both  the  exponential  and  Riehardson  regimes.  Boff¬ 
etta  and  Sokolov  (2002),  using  3D  direet  numerical  simula¬ 
tions  of  homogeneous  turbulenee,  emphasized  that 
deviations  from  the  Riehardson  regime  is  a  eonsequenee 
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of  turbulence  intermittency.  Balloon  experiments  in  the 
atmosphere  exhibited  an  exponential  regime  at  short  sepa¬ 
rations  (Morel  and  Lareheveque,  1974;  Er-El  and  Peskin, 
1981)  and  different  regimes  at  larger  scales  such  as  the  dif¬ 
fusive  law  (Morel  and  Lareheveque,  1974)  or  a  mixed 
regime  between  the  ballistic  and  Richardson  laws  (Er-El 
and  Peskin,  1981),  while  the  Richardson  regime  is  observed 
by  La co rata  et  al.  (2004)  for  separations  in  the  100— 
1000  km  range,  followed  by  a  diffusive  regime.  Experiments 
with  an  atmospheric  global  circulation  model  showed  all 
three  regimes:  exponential,  ballistic,  and  Richardson, 
depending  on  the  latitude  (Huber  ct  al,  2001). 

The  investigations  of  relative  dispersion  in  the  ocean  are 
relatively  few.  While  thousands  of  drifters  have  been 
released  during  the  WOCE  (World  Ocean  Circulation 
Experiment),  the  number  of  drifter  launches  that  are 
near-simultaneous  and  in  close  proximity  is  very  small. 
As  such,  this  large  data  set  is  not  suitable  for  a  systematic 
investigation  of  relative  dispersion  in  the  ocean.  Several 
process-oriented  ocean  experiments  have  been  conducted. 
Most  notably,  LaCasee  and  Bower  (2000)  estimated  rela¬ 
tive  dispersion  from  sub-surfaee  floats  released  during  sev¬ 
eral  observational  programs  in  the  North  Atlantic.  They 
found  a  diffusive  regime  in  the  eastern  part  of  the  domain 
and  a  Richardson  regime  in  the  western  part  of  the  domain 
that  was  attributed  to  an  inverse  turbulent  cascade  of 
energy.  The  exponential  regime,  as  expeeted  from  the  ens- 
trophy  cascade  and/or  ehaotie  adveetion,  was  not  resolved, 
possibly  beeause  the  initial  pair  separations  were  too  large. 
LaCasee  and  Ohlmann  (2003)  revisited  this  problem  using 
the  rich  drifter  data  set  in  the  Gulf  of  Mexico.  They  found 
an  exponential  regime  for  particle  separations  between  1 
and  4(1-50  km,  and  a  mixture  of  ballistic  and  Richardson 
regimes  up  to  scales  of  several  hundred  kilometers.  Laeora- 
ta  et  al.  (2001)  show  a  mixed  regime  between  ballistic  and 
Richardson  for  spatial  scales  of  40-200  km  using  drifters 
released  in  the  Adriatic  Sea.  Ollitrault  et  al.  (2005)  used 
another  set  of  sub-surfaee  floats  in  the  North  Atlantic 
and  documented  a  brief  exponential  regime  followed  by  a 
Riehardson’s  regime  for  separation  distances  between  40 
and  300  km  and  a  diffusive  regime  for  larger  separations. 

The  study  of  relative  dispersion  using  ocean  models  is 
surprisingly  limited  as  well.  Using  multi-layer,  quasi-gco- 
strophie,  double-gyre  simulations,  Berloff  et  al.  (2002) 
and  Berloff  and  McWilliams  (2002)  found  a  zonal  asymme¬ 
try  consistent  with  the  results  of  LaCasee  and  Bower 
(2000),  ranging  from  a  diffusive  regime  at  the  eastern  to  a 
ballistic  regime  at  the  western  boundaries.  Iudieone  et  al. 
(2002)  used  a  realistic,  eddy-permitting  (Ax  =  25  km)  circu¬ 
lation  model  of  the  Mediterranean  Sea  and  reported  an 
exponential  regime  at  meso-seales  and  a  ballistic  regime 
up  to  sub-basin  scales.  They  have  also  found  that,  while 
uncertainties  in  time  sampling  of  drifter  trajectories  lead 
to  differences  at  small  scales,  the  large-scale  dispersion 
remains  robust.  Similar  findings  are  reported  by  Braeco 
et  al.  (2004)  who  find  extremely  small  differences  in  relative 
dispersion  results  in  a  comparison  of  three-dimensional, 


quasi-geostrophie  turbulence  and  similar  slices  of  two- 
dimensional,  barotropic  flows. 

The  main  objective  of  the  present  study  is  to  contribute 
to  this  on-going,  long-standing  debate  on  the  classical 
problem  of  relative  dispersion,  particularly  given  the  rela¬ 
tively  small  number  of  studies  that  have  been  conducted 
directly  on  ocean  flows.  This  study  has  two  novel  aspects. 
First,  we  analyze  synthetic  drifter  trajectories  from  the 
Navy  Coastal  Ocean  Model  (NCOM)  with  high-resolution 
(1  km)  and  realistic  wind,  river  outflow,  tidal,  and  atmo¬ 
spheric  forcings.  The  model  has  been  shown  to  lead  to 
accurate  analyses  on  the  basis  of  operational  programs 
conducted  in  the  Adriatic  Sea  (Martin  et  al.,  2006;  Haza 
et  al,  2007).  The  Adriatic  circulation  provides  a  relatively 
extreme  ease  of  a  bounded  ocean  flow  dominated  by  coher¬ 
ent  features  including  intense  boundary  currents  and  sepa¬ 
ration  regions  and,  in  this  sense,  is  quite  far  removed  from 
a  homogeneous  or  isotropic  turbulence.  Second,  the  rela¬ 
tively  fine  spatial  and  temporal  model  resolution  allows 
us  to  address  questions  of  scale  dependence  and  velocity 
sampling  in  the  determination  of  model  Lagrangian  statis¬ 
tics.  As  discussed  by  Iudieone  et  al.  (2002),  there  is  an 
increasing  emphasis  on  Lagrangian  analyses  of  OGCM 
output.  Such  analyses  are  typically  conducted  off-line  with 
sub-sampled  time  slices  of  previously  computed  Eulerian 
model  velocity  fields.  Also,  remotely  sensed  altimetry  data 
are  often  analyzed  to  construct  approximate  velocity  fields 
for  the  purpose  of  Lagrangian  analyses.  The  observations 
necessarily  result  in  some  degree  of  both  spatial  and  tempo¬ 
ral  smoothing.  Questions  concerning  the  effects  of  temporal 
and  spatial  filtering  of  the  input  Eulerian  velocity  on  the 
resulting  Lagrangian  structures  and  statistics  can  be 
addressed  in  the  context  of  a  single  high-resolution  model 
sueh  as  the  one  considered  here. 

In  this  paper,  wc  address  three  main  questions: 

•  What  are  the  regimes  of  relative  dispersion  in  the  Adri¬ 
atic  Sea,  based  on  a  high-resolution  coastal  ocean 
model? 

•  How  arc  these  dispersion  regimes  modified  as  a  result  of 
uncertainties  in  the  Eulerian  velocity  field  as  represented 
by  their  spatial  and  temporal  smoothing? 

•  What  are  the  spatial  patterns  of  relative  dispersion  in  the 
model  Adriatic  and  how  robust  are  these  distributions  to 
perturbations  in  the  input  Eulerian  velocity  fields? 

It  is  also  worth  noting  that  spatial  averaging  affects  the 
turbulent  wave-number  spectrum,  while  time  averaging 
reduces  the  variability  of  the  Eulerian  field.  As  such,  a 
fourth  objective  of  this  study  is  to  make  progress  towards 
resolving  the  question  of  the  relative  importance  of  turbu¬ 
lent  cascade  processes  versus  chaotic  adveetion  arising 
from  the  time  variability  of  coherent  flow  features  on  the 
behavior  of  the  relative  dispersion. 

The  paper  is  organized  as  follows:  the  numerical  model 
is  described  in  Section  2,  the  Lagrangian  tools  are  defined 
in  Section  3,  results  for  the  full  field  and  filtered  velocities 
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arc  presented  in  Section  4,  and  the  main  results  are  summa¬ 
rized  in  Section  5. 

2.  Coastal  model  and  configuration 

The  ocean  model  used  for  this  study  is  the  NCOM  as 
described  in  Martin  (2000),  with  some  improvements  as 
described  in  Morey  et  al.  (2003)  and  Barron  et  al.  (2006). 
This  is  a  hydrostatic  model,  which  is  similar  in  its  physics 
and  numerics  to  the  Princeton  Ocean  Model  (POM) 
(Blumberg  and  Mellor,  1987),  but  uses  an  implicit  treat¬ 
ment  of  the  free  surface  and  a  hybrid  vertical  grid  with 
sigma  coordinates  in  the  upper  layers  and  (optionally)  level 
coordinates  below'  a  user-speeified  depth. 

The  model  equations  include  a  source  term  that  can  be 
used  for  river  inflows.  A  third-order  upwind  method  (Hol¬ 
land  et  al ,  1998)  was  used  for  advection.  Vertical  mixing 
was  computed  using  the  Mellor-Yamada  Level  2  seheme 
(Mellor  and  Yamada,  1974).  The  equation  of  state  used 
was  that  of  Mellor  ( 1991 ). 

The  ocean  model  domain  consists  of  the  entire  Adriatic 
Sea,  a  sub-basin  of  the  Mediterranean  Sea,  and  it  includes 
the  Strait  of  Otranto  and  a  small  part  of  the  northern 
Ionian  Sea.  The  setting  of  the  model  for  the  Adriatic  Sea 
is  described  in  detail  in  Martin  et  al.  (2006).  The  horizontal 
grid  resolution  is  1019.5  m.  The  vertical  grid  consists  of  32 
total  layers,  with  22  sigma  layers  used  from  the  surface 
down  to  a  depth  of  291  m  and  level  coordinates  used  below 
291  m.  Hence,  the  grid  is  like  a  regular  sigma  coordinate 
grid  in  water  shallower  than  291  m  and  similar  to  a  level 
grid  in  deeper  water.  The  vertical  grid  is  uniformly 
stretched  from  the  surface  downward  with  a  maximum 
thiekness  of  the  upper  layer  of  2  m  and  a  maximum  depth 
of  1262  m 

Initial  conditions  and  daily  boundary  conditions  were 
taken  from  a  hind-cast  of  a  global  model  (Barron  et  al., 
2004).  The  numerical  treatment  of  the  boundary  conditions 
includes  the  Flather  radiation  condition  (Flather  and  Proc¬ 
tor,  1983)  for  the  surfaee  elevation  and  depth-averaged 
normal  velocity,  Orlanski  radiation  conditions  (Orlanski, 
1976)  for  the  tangential  velocities  and  sealar  fields,  and  a 
relaxation  to  the  temperature  and  salinity  from  the  global 
model  near  the  open  boundary.  Tidal  forcing  was  provided 
using  tidal  elevation  and  depth-averaged  normal  and  tan¬ 
gential  velocities  at  the  open  boundaries  from  the  Oregon 
State  University  (OSU)  tidal  data  bases,  which  are  derived 
from  satellite  altimetry  data  (Egbert  and  Erofeeva,  2003). 
Data  from  the  OSU  Mediterranean  tidal  data  base  were 
used  for  the  Kl,  Ol,  M2,  and  S2  constituents  and  data 
from  the  OSU  global  data  base  were  used  for  PI,  Ql, 
K2,  and  N2.  Tidal  potential  forcing  for  these  eight  constit¬ 
uents  was  used  in  the  interior  of  the  model  domain. 

Atmospheric  forcing  was  obtained  from  the  Coupled 
Ocean/Atmosphere  Meso-scale  Prediction  System 
(COAMPS)  (Hodur,  1997).  The  COAMPS  setup  for  the 
Adriatic  consists  of  a  triply  nested  grid  with  resolutions 
of  36,  12,  and  4  km  (Pullen  et  al.,  2003).  The  outer  grid 


of  this  nested  grid  system  eovers  most  of  Europe  and  the 
Mediterranean  and  the  inner  4-km  grid  eovers  the  entire 
Adriatic  and  part  of  the  Tyrhennian  Sea.  COAMPS  itself 
is  nested  within  the  Navy  Operational  Global  Atmospherie 
Prediction  System  (NOGAPS)  (Rosmond  et  al.,  2002). 
River  and  runoff  inflows  for  the  Adriatic  were  taken  from 
the  monthly  elimatological  data  base  of  Raicich  (1996). 
This  data  base  includes  discharges  for  about  39  rivers 
and  runoff  inflows  along  a  number  of  sections  of  the  Adri¬ 
atic  coastline.  Raicich’s  monthly  climate  values  were  used 
for  all  the  inflows,  exeept  that  daily  observed  discharge  val¬ 
ues  were  used  for  the  Po  River  (Rich  Signell,  personal 
communication). 

This  model  configuration  has  shown  good  agreement 
with  tidal  and  de-tided  ADCP  current  observations  and 
reproduced  some  of  the  observed  temperature  and  salinity 
spatial  profiles  (Martin  et  a.l„  2006).  Additionally,  the 
model  response  to  Bora  wind  events  displayed  the  usually 
observed  pattern  of  cyclonic  and  anticyclonic  gyres. 

This  model  is  appropriate  for  the  present  study  for  two 
main  reasons.  First,  a  recent  operational  experiment  aimed 
at  targeted  release  of  drifters  in  high  dispersion  regions  in 
the  Adriatic  Sea  using  the  same  model  configuration 
showed  the  ability  of  NCOM  to  reproduce  many  aspeets 
of  the  meso-scale  circulation  in  real  time  (Haza  et  al., 
2007).  Second,  the  high  spatial  and  temporal  resolution 
of  the  model  combined  with  a  rich  variety  of  forcing  mech¬ 
anisms  and  internal  coastal  flow  dynamics  lead  to  a  surface 
velocity  field  exhibiting  features  at  a  wide  range  of  seales 
and  with  high  temporal  variability. 

A  snapshot  of  the  velocity  field  in  the  full  numerical 
domain  is  depicted  in  Fig.  la.  The  reader  is  referred  to  Pou- 
lain  (1999,  2001),  Cushman-Roisin  ct  al.  (2001,  2007),  Fal- 
co  et  al.  (2000),  Maurizi  et  al.  (2004)  and  Martin  et  al. 
(2006)  for  a  comprehensive  diseussion  of  the  circulation 
features  and  data  sets.  The  main  point  is  that  the  NCOM 
simulation  employed  here  adequately  approximates  the 
main  features  of  the  upper  circulation  in  the  Adriatic  Sea. 

3.  Experimental  setup 

The  bulk  of  the  results  that  follow  concern  the  disper¬ 
sion  statistics  of  synthetic  drifters  advected  for  30-60  days 
using  the  NCOM  surface  velocity  field  of  September  and 
October  2002.  Throughout,  we  use  a  standard,  fourth- 
order,  Runge-Kutta  seheme  to  solve  the  Lagrangian 
equations  with  an  integration  time-step  of  At  =  1  h  corre¬ 
sponding  to  the  model  archive  time.  Spatial  interpolation 
is  achieved  by  a  third-order  polynomial  routine. 

A  subset  of  partiele  trajectories  is  displayed  in  Fig.  lb; 
they  describe  the  cyclonic  circulation  of  the  Adriatic  Sea 
corresponding  to  the  Western  and  Eastern  Adriatic  Cur¬ 
rents  (WAC  and  EAC,  respectively),  including  three  other 
cyclonic  gyres  within  the  main  cyclonic  circulation.  The 
velocities  of  the  boundary  currents  (in  particular  the 
WAC)  are  of  the  order  of  50-100  cm/s  (as  seen  in  the  snap¬ 
shot  of  the  surfaee  velocity  field  in  Fig.  la),  while  those  of 
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Fig.  1 .  (a)  Snapshot  of  the  raw  NCOM  surface  velocity  field,  (b)  One-month  long  trajectories  of  a  subset  of  74  particles  advccted  with  the  NCOM  velocity 
field  corresponding  to  September  2002. 


the  Adriatic  interior  are  ~  10-20  cm/s  depending  on  the 
surface  winds.  In  the  regions  away  from  the  boundary  cur¬ 
rents,  the  trajectories  display  loops  with  radii  of  1-2  km, 
indicating  the  presence  of  inertial  oscillations.  Perhaps 
the  most  striking  feature  of  the  trajectories  from  Pig.  lb 
is  the  degree  of  control  exerted  by  the  strong  boundary  cur¬ 
rents,  with  large  transport  in  the  longitudinal  direction  and 
relatively  limited  cross-basin  exchange. 

Given  the  need  to  adequately  sample  the  large  degree 
of  spatial  inhomogencity  in  the  Adriatic  model  and  the 
desire  to  accurately  map  the  scale-dependence  of  the  drif¬ 
ter  dispersion  across  the  complete  range  of  available 
model  length  scales,  we  consider  two  distinct  initial 
launch  grids  and  their  corresponding  statistical  sampling 
strategies.  First,  a  relatively  dense  array  of  drifters  is  ini¬ 
tialized  on  a  uniform  regular  grid  with  5.1  km  spacing. 
This  corresponds  to  a  total  of  N  =  5172  drifters,  placing 
an  upper  bound  of  N(N  —  l)/2  ~  13  x  106  possible  drifter 
pairs  from  which  to  compute  statistics.  The  second  config¬ 
uration  is  composed  of  M  =  5N  drifters  initialized  in  N 
crosses  with  separation  Aa  ~  1  km  centered  on  the  origi¬ 


nal  N  launch  locations.  In  this  case,  only  the  statistics 
of  the  fixed  number  of  particle  pairs  with  initial  separa¬ 
tion  Z)0  =  A  a  are  considered. 

As  mentioned  in  Section  1,  both  temporal  and  spatial  fil¬ 
terings  are  applied  to  the  Eulerian  velocity  fields,  with  the 
objective  of  not  only  assessing  the  importance  of  uncertain¬ 
ties  in  the  input  velocities  on  the  relative  dispersion,  but 
also  for  examining  the  role  of  turbulent  cascade  processes 
versus  chaotic  advection.  The  temporal  filtering  of  the  flow 
is  conducted  by  locally  averaging  the  velocity  field  over 
time-windows  of  varying  duration.  For  a  time  window 
given  by  T  =  MAt ,  the  filtered  velocity  u}  '  at  location  ij 
and  time  t  =  kAt  is 

,  m  k+M/2 

«fj(*)=77  uiAm)-  (5) 

M  m  k  M/2 

The  spatial  filtering  of  the  flow  is  carried  out  via  a  convo¬ 
lution  product  of  the  velocity  field  with  a  Gaussian  func¬ 
tion.  The  spatially  filtered  velocity  wFs  at  location  iQJ0 
and  time  t  is 
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«*F,0Wo)  = 


!(<-'(>)'  C/  /o ) 2 1 /« ' 


.eW-W+U-JtfV* 


(6) 


This  corresponds  to  a  Gaussian-weighted,  spatially  aver¬ 
aged  smoothing  of  the  horizontal  gradients  of  the  flow  with 
a  length  scale  ~2cx,  while  preserving  the  time  dependence. 
The  convolution  (6)  is  earned  out  with  a  —  0.5;  1 ;  2.5;  5; 
10;  20  (approximately  in  kilometer,  sinee  the  horizontal 
grid  spaeing  is  1  km). 


4.  Results 

4.1.  Dispersion  from  the  unfiltered  veloc  ity  field 

Relative  dispersion  is  first  calculated  from  the  trajectory 
ensemble  obtained  from  the  raw  (unfiltered)  velocity  field 
using  the  complete  set  of  5172  synthetic  drifters  and 
all  possible  dnfter  pair  combinations.  The  anisotropy  of 
dispersion  is  well  known  in  the  Adriatic  Sea  (Falco  et  al , 
2000;  Laeorata  et  al.,  2001).  As  sueh,  the  relative  dispersion 
is  calculated  in  the  zonal  and  meridional  directions  sepa¬ 
rately.  Fig.  2  shows  the  relative  dispersion  curves  based 
on  trajectories  with  the  original  launch  spaeing  of 
D0  =  5.1  km.  The  relative  dispersion  eurves  are  fairly  well 
behaved  due  to  the  large  number  of  synthetic  drifter  pairs 
included  in  the  ealeulation,  amounting  to  more  than  9800 
pairs  initially,  and  down  to  about  6500  pairs  at  the  end 
of  the  2-month  experiment.  For  times  longer  than  a  few 
days,  relative  dispersion  in  the  meridional  direction  is  best 
approximated  by  a  t7  growth,  whereas  relative  dispersion  in 
the  zonal  direction  shows  significantly  slower  growth  seal¬ 
ing  approximately  like  z1,3.  The  total  relative  dispersion 
scales  as  D2{t)  ~  tl  9.  Sueh  a  super-diffusive  or  nearly  ballis¬ 
tic  regime  is  typically  associated  with  shear  dispersion  (Iud- 
icone  et  al.,  2002;  LaCasee  and  Ohlmann,  2003)  and 
persistent  velocity  correlations.  Both  effects  arc  clearly  evi¬ 
denced  in  the  Adriatic  Sea  eireulation  due  to  the  presence 
of  boundary'  currents  (Fig.  1).  As  the  direction  of  the 
boundary  currents  is  imposed  by  the  morphology  of  the 
Adriatic  coastline,  the  similar  power  laws  of  D2(t)  and 
indicate  that  the  relative  dispersion  is  enhanced  along 
the  main  axis  of  the  Adriatic  basin.  For  times  longer  than  a 
month,  relative  dispersion  in  both  directions  slows  down 
considerably  to  a  sub-diffusive  regime  of  t°  .  This  regime 
corresponds  to  dispersion  at  scales  comparable  to  the  basin 
size  (Artalc  et  al.,  1997). 

Relative  dispersion  as  a  function  of  separation  distance 
ean  be  quantified  naturally  using  the  FSLE.  While  it  is 
generally  true  that  the  FSLE  defined  by  Eq.  (4)  is  insen¬ 
sitive  to  the  choice  of  the  parameter  y.  numerical  and  sta¬ 
tistical  considerations  dietate  the  range  of  this  number  in 
applications.  The  minimum  value  of  y  is  set  by  the  tempo¬ 
ral  resolution  At  of  the  trajectories  sueh  that  the  time 
required  for  partiele  pairs  to  disperse  from  <50  to  a<50  be 
longer  than  At.  Thus,  a<50  -  30  >  AtAi\  where  Ar  is  the 
velocity  difference  between  partieles  in  a  pair.  This  leads 


Fig.  2.  (a)  Relative  dispersion  />2(/)  of  the  raw  field  computed  in  the  zonal 
and  meridional  directions  from  the  synthetic  trajectories  with  original 
initial  release  distance  of  Do  5. 1  km.  (b)  Sensitivity  of  FSLE,  to  the 
time  step  A /  and  parameter  a  for  the  same  set  of  synthetic  drifter 
trajectories. 


to  am«n  —  I  +  AtAv/S0.  In  this  experiment.  At  =  1  h, 

=  1  km,  and  an  estimate  of  Av  ~  3  km/day  given  by 
the  Lagrangian  structure  function  computed  from  the 
synthetic  trajectories  at  the  minimum  seale  of  <5=1  km. 
Therefore,  amin~1.16  is  the  estimated  minimum  value 
for  accurate  ealeulation  of  the  FSLE. 

Fig.  2b  displays  plots  of  A(<5)  for  various  values  of  a  in 
the  range  of  1.05  1.3.  The  FSLE  eurves  collapse 

for  a  >  1.2,  in  agreement  with  the  above  estimate.  There 
is  an  increasing  flattening  of  the  FSLE  eurves  at  seales 
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smaller  than  10  km  for  a  <  1.13,  ultimately  leading  to  a 
spurious  exponential  regime.  The  length  of  the  exponential 
regime  for  x  <  amin  is  further  enhanced  if  the  velocity  field 
is  sub-sampled  every  3  h.  For  the  following  experiments, 
we  assign  x  =  1  2,  since  it  is  preferable  to  use  a  value  of  x 
as  small  as  possible  to  resolve  the  scale  separation  (Laeora- 
ta  et  al.,  2001),  and  this  is  the  smallest  x  in  this  set  of  exper¬ 
iments  where  k  is  insensitive  to  the  precise  value  of  the 
parameter. 

Focusing  on  the  nearly  identical  curves  of  A(<5)  ^  1.2  in 
Fig.  2,  we  note  the  following.  For  scales  larger  than  20  km, 
the  FSLE  follows  a  power  law  /(<$)  ~  S  °'8.  This  super-dif¬ 
fusive  regime  is  fairly  consistent  with  D2(t)  ~  t 2  "  (which 
corresponds  to  k  ~  S  °'9)  seen  in  Fig.  2a.  The  two  metrics, 
D2(t)  and  k(S)  do  not  necessarily  lead  to  identical  results 
due  to  sampling  differences  (Laeorata  et  al.,  2001;  LaCasee 
and  Ohlmann,  2003).  The  relative  dispersion  curve  aver¬ 
ages  typically  have  increasing  sample  sizes  for  smaller  sep¬ 
arations,  while  the  FSLE  curves  result  from  separation 
times  averaged  over  larger  numbers  of  particle  pairs  at  lar¬ 
ger  scales.  As  such,  the  small  discrepancy  in  sealing  seen 
here  appears  satisfactory.  The  sealing  relation  for  the 
model  FSLE,  /  ~  S  0  8  for  <5  ^  20  km,  is  in  very  good 
agreement  with  k  %  S  0  84  in  the  range  of  30  <  S  <  200  km 
obtained  by  Laeorata  et  al.  (2001,  their  Fig.  4)  on  the  basis 
of  trajectory  data  from  37  drifters  launched  in  the  Adriatic, 
although  the  values  of  k  for  the  synthetic  trajectories  are 
higher. 

The  Adriatic  circulation  is  known  to  exhibit  strong  sea¬ 
sonality  due,  in  part,  to  changes  in  the  form  of  the  prevail¬ 
ing  wind  stress.  In  order  to  assess  the  effect  of  seasonality 
on  the  relative  dispersion  and  determine  how  representative 
any  single-month  analysis  may  be,  we  examine  8  months  of 
available  NCOM  output  This  includes  three  months  in 
2002  (non-assimilated)  and  five  months  of  sea-surfaee  tem¬ 
perature  assimilating  results  from  2006.  Fig.  3  displays 
1 -month  records  of  D2(t)  and  A(<5)  for  each  of  these  8  months. 
Note  that  in  spite  of  the  summer-fall  predominance  in  the 
data  set,  the  relative  dispersion  at  all  times  and  scales  for 
March  2006  (winter-spring)  and  June  2006  (spring-sum¬ 
mer)  does  not  differ  significantly  from  the  other  months. 
We  obtain  power  laws  of  D2(i)  varying  between  f1  9  and 
/2  3  (with  an  average  ~/2'2)  for  D0  =  1  km,  and  of  /(<$) 
between  S  0  *  and  S  09  (with  an  average  of  °  8)  for 
S  >  20  km.  As  seen  in  the  FSLE  curves,  there  are  small 
but  distinguishable  differences  between  the  data  assimilat¬ 
ing  and  non-data  assimilating  fields  and  a  marked  increase 
in  the  monthly  variability  within  the  data  assimilating  set. 
The  relatively  robust  sealing  indicates,  however,  that  the 
dynamics  which  control  the  relative  dispersion  are  present 
in  all  seasons,  and  in  both  configurations  of  the  model. 

Perhaps  the  most  interesting  aspect  of  any  of  the  FSLE 
curves  computed  from  the  uniform  D0  =  5.1km  drifter 
array  using  all  possible  particle  pair  statistics  is  the  lack 
of  a  clear  exponential  regime  (constant  /(£))  at  small  sepa¬ 
ration  scales.  Such  a  plateau  would  be  expected  in  any  eha- 
otie  adveetion  regime  where  the  dispersion  is  dominated  by 


Fig.  3.  (a)  Relative  dispersion  D2{t)  for  f)0  I  km  (chance  pairs)  and  (b) 
FSLE  /(<$)  curves  for  all  available  monihs  in  2002  and  2006.  The  range  in 
power  laws  is  displayed  by  the  black  lines,  while  the  average  power  law  is 
indiealed  in  red. 


the  time  dependence  of  a  spatially  smooth  velocity  field.  In 
the  present  context,  however,  there  seems  to  be  no  clear  law 
governing  dispersion  at  scales  in  the  range  of 
1  <  S  <  20  km  The  reasons  for  this  almost  surely  involve 
the  marked  inhomogeneity  in  the  Adriatic  circulation  and 
the  dominating  presence  of  the  strong  western  boundary 
current.  Strong  differences  in  the  statistics  of  ‘interior’ 
and  ‘boundary  current’  drifter  pairs  are  coupled  with  the 
highly  non-uniform  rates  at  which  these  distinct  regions 
are  sampled  by  any  initial  distribution  of  drifters. 

The  effects  of  this  are  clearly  seen  in  Fig.  4  where  the  dif¬ 
ferences  between  the  D0  =  5.1  km  launch  using  all  possible 
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Fig.  4.  (a)  Total  relative  dispersion  DV)  with  0{)  I  km  computed  from  chance  pairs  (black  crosses)  and  from  the  original  pairs  of  the  ‘pair-launch’ 
selling  (gray  crosses),  (b)  Corresponding  FSLE  curves  for  the  two  sampling  strategies.  Inset  in  upper  figure  shows  the  location  of  all  chance  pairs  wilh 
I  km  during  lhe  observation  period. 


particle  pairs  and  the  D0=  1.0  km  launch  using  only  the 
original  pair  statistics  are  compared  for  September  2002. 
In  general,  it  is  possible  to  calculate  relative  dispersion 
based  not  only  on  the  particles  that  are  launched  together 
(so-called  original  pairs),  but  also  based  on  those  (so-called 
chance  pairs)  that  approach  one  another  at  some  later  time 
(LaCasce  and  Bower,  2000).  Chance  pair  sampling  is  typi¬ 
cally  employed  to  increase  the  sample  size  when  only  a  lim¬ 
ited  number  of  drifter  trajectories  are  available.  This  is 
usually  the  ease  in  oceanic  observations,  but  not  a  critical 
factor  when  synthetic  drifter  data  can  be  generated  from 
numerical  model  output.  Nevertheless,  using  chance  pairs. 


it  is  possible  to  calculate  the  time  evolution  of  the  relative 
dispersion  of  particles  starting  from  smaller  initial  dis¬ 
tances  than  those  of  the  original  pairs,  namely  down  to 
the  numerical  model  grid  size  of  1  km. 

Fig.  4a  shows  that  relative  dispersion  for  the  7)0  =  1  km 
pair-launch  displays  D2{t)  ~  rA  which  is  somewhat  steeper 
than  that  of  the  ehanee  pair  computation  with  larger  initial 
separation,  although  the  statistics  at  long  times  are  consid¬ 
erably  less  robust  in  this  ease.  More  striking  differences 
between  the  two  sampling  strategies  are  seen  in  the  FSLE 
results  shown  in  Fig.  4b.  While  the  behavior  is  unchanged 
at  scales  larger  than  20  km,  the  fixed  sampling  pair  strategy 
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exhibits  a  distinct  exponential  regime  at  smaller  scales.  The 
reasons  for  this  difference  are  apparent  when  one  examines 
the  spatial  locations  of  all  chance  pairs  approaching  within 
Z)0  =  1  km  as  shown  in  the  inset  of  Fig.  4a,  The  inhomoge¬ 
neity  of  the  flow  field  and  the  lack  of  any  ergodieity  in  the 
drifter  trajectories  implies  that  chanee  pairs  with  small  sep¬ 
arations  occur  with  much  higher  frequency  on  the  energetie 
boundaries  of  the  Western  Adriatic  Current.  At  small 
scales,  the  chanee  pair  statistics  preferentially  weight  this 
flow  feature.  In  contrast,  fixing  the  number  of  particle  pairs 
at  the  outset,  as  in  the  original  pair  strategy,  preferentially 
weights  less  energetic  interior  gyre  flow  features  which  pro¬ 
duce  larger  separation  times. 

The  controlling  effeet  of  energetic  WAC  also  explains 
the  marked  sensitivity  of  the  computed  FSLE  curves  to 
the  choice  of  doubling  parameter  y.  Choosing  y  too  small 
effectively  removes  those  particle  pairs  in  the  WAC  which 
separate  the  smallest  distances  on  time  scales  faster  than 
the  At  imposed  by  the  model  archiving.  As  expected,  aver¬ 
aged  velocity  differences  computed  from  the  Lagrangian 
structure  function  for  the  fixed  pair  launch  are  consider¬ 
ably  less  (Ai?~  1  km/day  at  1  km)  than  those  computed 
from  the  chance  pairs  and  the  corresponding  /(<$)  can  be 
reliably  computed  with  even  smaller  values  of  the  doubling 
parameter,  amm. 

The  results  in  Fig.  4  indicate  that  the  relative  dispersion 
in  the  model  Adriatic  on  scales  <20  km  is  extremely  sensi¬ 
tive  to  both  the  model  resolution  of  the  WAC  and  to  how 
well  this  region  is  sampled  by  the  Lagrangian  drifters.  The 
dynamics  of  the  synthetic  drifters  at  such  scales  are,  in  turn, 
dependent  upon  the  details  of  the  spatial  and  temporal 
interpolation  schemes  used.  All  processes  at  these  scales 
are  likely  to  be  sensitive  to  the  spatial  and/or  temporal 
accuracy  of  the  model.  Therefore,  it  is  at  these  scales  where 
one  would  anticipate  the  effeet  of  the  uncertainties  by  spa¬ 
tial  and  temporal  smoothing  to  be  most  significant.  In  what 
follows,  we  restrict  ourselves  to  the  chanee  pair  launch 
strategy  which  makes  use  of  all  possible  drifter  pair  statis¬ 
tics  from  uniform  launches  of  N  =  5172  synthetic  drifters. 
This  choice  provides  large  numbers  of  sample  pairs  at  all 
scales  while  mimicking  the  observational  reality  of  allowing 
the  Lagrangian  dynamics  to  control  the  location  of  individ¬ 
ual  drifter  pairs  at  any  time. 

4.2.  Effect  of  temporal  averaging  on  relative  dispersion 

Having  quantified  relative  dispersion  and  FSLE  from 
the  raw  field,  we  now  look  at  the  effect  of  temporal  aver¬ 
aging,  First,  snapshots  of  the  flow  field  and  synthetic  tra¬ 
jectories  are  plotted  in  Fig.  5  in  order  to  assess  how 
averaging  the  Eulerian  field  impacts  the  drifter  trajectories. 
Three  levels  of  temporal  filtering  are  employed:  1,5,  3,  and 
30  days.  A  1.5-day  filtering  is  commonly  used  for  the 
removal  of  inertial  oscillations  exhibited  by  the  drifter  tra¬ 
jectories  in  the  Adriatic  Sea.  The  subset  of  trajectories 
after  sueh  filtering  of  the  flow-field  indeed  shows  that  most 
of  the  loopers  present  in  the  Adriatic  interior  (Fig.  lb) 


have  been  eliminated  (lower  panel.  Fig.  5b).  We  note  that 
the  3-day  moving  average  reduces  by  half  the  intensity  of 
the  anticyclonic  eddy  in  the  southern  Adriatic,  due  to  its 
propagation,  while  the  overall  structure  of  the  boundary 
currents  are  fairly  well  preserved  in  the  three  smoothing 
cases.  Eddies  disappear  almost  entirely  from  the  flow  field 
in  the  ease  of  30-day  temporal  filtering.  Similar  consider¬ 
ations  apply  to  other  fast-moving  events,  while  the 
quasi-permanent  coherent  features  of  the  flow  field  remain 
after  time  averaging.  Such  a  flow  field  leads  to  extremely 
smooth  trajectories  that  are  mainly  confined  to  the  bound¬ 
ary  currents  and  exhibit  significantly  less  zonal  (cross- 
basin)  transports.  This  implies  that  much  of  the  zonal 
transport  in  this  Adriatic  Sea  simulation  is  carried  out 
by  episodic  events,  sueh  as  meandering  of  the  boundary 
currents,  swirls,  filaments,  and  eddies.  As  shown  in  Seetion 
4.4,  the  meandering  and  eddying  of  the  WBC  is  a  major 
source  of  rapid,  meso-seale  and  submeso-scale  partiele  dis¬ 
persion  in  the  Adriatic  basin. 

The  relative  dispersion  and  FSLE  curves  for  velocity 
fields  subjected  to  temporal  averaging  are  plotted  in 
Fig.  6a  and  b,  respectively.  The  eurves  from  the  raw  fields 
are  also  shown  for  comparison.  From  both  /)“(/)  and  A((5), 
it  is  apparent  that  time-filtering  influences  the  small-scale 
dispersion  only.  In  the  case  of  Z)2(/),  the  influence  of  the 
removal  of  time  variability  manifests  itself  by  an  increased 
delay  for  the  onset  of  dispersion  in  the  range  of  a  few  hours 
to  2  days.  The  small-scale,  short-time  dispersion  is  reduced. 
The  large-scale  dispersion  remains  unchanged  from  the 
super-diffusive  regime  of  D2{t)  ~  tl  9  In  the  case  of  the 
FSLE,  the  effeet  of  time  smoothing  manifests  itself  by 
reduced  dispersion  (smaller  FSLE)  at  scales  in  the  range 
of  1-20  km,  which  essentially  covers  the  submeso-scale 
and  meso-seale  variability  that  is  filtered  out.  Dispersion 
at  the  smallest  scale  computed  here,  8  =  1  km,  is  reduced 
by  half  using  the  monthly  mean  velocity  field.  As  in 
Z)2(/),  dispersion  at  larger  seales  remain  unchanged,  dis¬ 
playing  /  ~  8  0  8  for  8  >  20  km. 

We  note  that  the  features  of  the  circulation  remaining  in 
the  mean  flow  are  the  boundary  currents  and  basin-scale 
gyres  after  30-day  temporal  smoothing  (upper  panel, 
Fig.  5c).  Wc  conclude  that  the  mean  boundary  currents 
and  persistent  stretching  at  the  interface  between  the  basin 
gyres  are  actually  the  main  circulation  features  controlling 
the  dispersion  at  large  scales.  These  features  are  not 
affected  by  the  temporal  filtering  of  the  flow  field.  As  such, 
the  D2(t)  ~  tl  9  or  /.(£)  ~  8 “°  8  regimes  appear  to  be  a  con¬ 
sequence  of  horizontal  shear  dispersion  and  persistent  iso¬ 
lated  regions  of  high  strain  in  the  Eulerian  field. 

4.3.  Effect  of  spatial  smoothing  on  relative  dispersion 

The  effect  of  spatial  smoothing  on  the  relative  dispersion 
is  investigated  next  by  applying  a  Gaussian  filter  to  the 
velocity  field  for  different  values  of  a.  The  upper  panel  of 
Fig.  7  illustrates  the  progressive  smearing  of  the  circulation 
features  with  a  values  of  2.5,  5,  and  10  km  The  internal 
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Fig  5.  Upper  panel,  snapshots  of  the  surface  velocity  field  after  applying  time  smoothing  with  a  window  size  of  (a)  1.5  days,  (b)  3  days,  and  (c)  30  days 
The  colors  indicate  the  speed  (in  cm/s).  Lower  panel,  trajectories  of  74  synthetic  drifters  advecled  with  the  same  low-pass  filtered  velocity  fields. 


deformation  radius  in  the  Adriatic  can  be  as  short  as  5  km 
(Cushman-Roisin  ct  al.,  2007).  As  such,  mcso-scalc  features 
like  the  anticyclomc  eddy  downstream  of  the  Gargano 
Cape  arc  smoothed  out  with  a  =  10  km  and  the  boundary 
currents  become  wider  and  slower  for  larger  values  of  cj , 


The  corresponding  trajectories  (lower  panel.  Fig.  7) 
show  how  different  the  effect  of  spatial  filtering  is  from  that 
of  temporal  filtering.  The  difference  in  the  trajectories  from 
the  spatially  filtered  velocity  fields  and  those  from  the  raw 
field  (Fig.  lb)  is  barely  apparent;  inertial  oscillations  arc 
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Fig.  6,  (a)  Relative  dispersion  D2{t)  and  (b)  FSLE  A(<5)  curves  in  cases  of 
lemporal  low-pass  fillcring  of  the  velocity  field. 


still  present  in  the  regions  of  low  kinetic  energy,  even  for 
large  <r,  and  cross-basin  transports  remain  virtually 
unmodified.  The  persistence  of  inertial  oscillations  is  not 
surprising  as  their  dynamical  balance  requires  only  time 
dependence  of  the  underlying  field  and  the  effect  of 
rotation 

D2  and  /  computed  for  six  different  values  of  a  ranging 
from  0.5  to  20  km  are  plotted  in  Fig.  8.  For  these  eases,  the 
trajectories  of  particles  venturing  too  close  to  the  coast 
were  discarded  to  avoid  the  spurious  dispersion  generated 
by  the  spatial  smearing  of  the  surface  boundary  layer  with 
adjacent  model  land  points.  A  condition  on  a  depth  shal¬ 
lower  than  20  m  was  sufficient  for  most  of  the  Adriatic 


basin  with  the  exception  of  the  Croatian  coast  where  the 
entire  circulation  around  the  islands  was  dismissed.  The 
FSLE  and  relative  dispersion  of  the  full  field  were  also 
recalculated  under  the  same  boundary  conditions  for 
consistency. 

For  0.5  ^  (T  ^  5  km,  the  spatial  filtering  aets  on  the  rel¬ 
ative  dispersion  in  a  similar  way  as  the  temporal  filtering, 
by  affecting  mostly  the  scales  below  10-20  km.  In  particu¬ 
lar,  (T  =  2.5  km  and  1 .5-day  low-pass  filtering  yield  very 
similar  curves.  We  conclude  that  spatial  smoothing  in  this 
range  eliminates  fine-scale  turbulent  features  that  are  also 
fast  enough  to  be  affected  by  the  1.5-day  temporal  filtering. 
As  shown  in  Section  4.4,  such  small  scale  structures  are 
predominantly  associated  with  instabilities  in  the  WAC 
and  eddies  within  the  Croatian  island  chain.  Note  also  that 
the  D 1  power  law  remains  the  same  in  this  range  of  a. 

For  c  ^  10  km,  the  large-scale  dispersion  starts  to  be 
affected,  which  corresponds  to  the  smearing  of  the  WAC 
and  EAC  observed  in  Fig.  7c.  The  dispersion  curves  D2(t) 
for  the  large  <r,  as  seen  in  Fig.  8a,  no  longer  display  a  power 
law.  The  FSLE  /  exhibits  now  an  exponential-type  regime 
for  5  ^  6  ^  20  km.  Therefore,  it  appears  that  spatial 
smoothing  large  enough  to  weaken  the  features  of  the  mean 
flow  influences  relative  dispersion  at  all  scales.  Relative  dis¬ 
persion  tends  to  be  dominated  by  chaotic  advection  at 
small  scales  since  time  variability  is  preserved  and  large- 
scale  relative  dispersion  tends  to  decrease  because  of 
reduced  horizontal  shears  in  the  boundary  currents  and 
wcakening/climination  of  high  rate  of  strain  regions 
throughout  the  domain. 

4.4.  Analysis 

One  of  the  striking  results  of  these  calculations  is  the 
robustness  of  the  large-scale  relative  dispersion  to  temporal 
and  spatial  filtering  of  the  velocity  field,  or  in  other  words, 
the  degree  of  control  exerted  by  the  large-scale  circulation 
on  the  traecr  evolution.  We  now  provide  arguments  to  help 
explain  this  finding  based  on  both  classical,  namely  energy 
cascade,  and  recent  dynamical-system  methods. 

Bennett  (1984)  put  forth  the  concept  of  two  dynamical 
regimes:  a  loeal  regime  where  the  tracer  evolution  at  a 
given  length  scale  is  controlled  by  the  velocity  field  at  the 
same  length  scale  and  a  non-local  regime  where  the  tracer 
evolution  is  controlled  by  the  velocity  field  at  much  larger 
scales.  In  particular,  Bennett  (1984)  showed  that  the  transi¬ 
tion  between  the  local  and  non-loeal  regimes  depends  on 
the  slope  of  the  kinetic  energy  spectrum.  When  the  slope 
is  steeper  than  —3,  then  non-loeal  control  of  the  tracer  field 
becomes  dominant.  While  this  analysis  relies  on  the  con¬ 
nection  between  the  velocity  structure  function  and  the 
energy  spectrum  for  homogeneous,  isotropic  turbulence, 
investigations  by  Babiano  et  al.  (1985)  and  Shepherd 
et  al.  (2000)  for  derived  atmospheric  velocity  data  further 
support  the  result.  While  the  Adriatic  model  produces  a 
velocity  field  is  distinctly  inhomogeneous  and  non-isotro- 
pie,  computations  of  the  wave-number  speetra  (not  shown) 
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1  ig.  7  Upper  panel,  snapshots  of  the  surface  velocity  field  after  applying  a  Gaussian  spatial  filter  with  (a)  a  —  2.5  km,  (b)  a  5  km,  and  (c)  a  10  km. 
Lower  panel,  trajectories  of  74  synthetic  drifters  advected  with  the  same  spatially  filtered  velocity  fields. 


indicate  spectral  slopes  significantly  steeper  than  k  \  The 
observed  ‘non-local1  behavior  of  the  relative  dispersion, 
dominated  by  the  presence  of  highly  sheared  boundary  cur¬ 
rents  and  large-seale  gyres,  is  at  least  not  inconsistent  with 
spectral  predictions. 


Next,  we  investigate  the  spatial  distribution  of  the 
FSLEs  and  how  temporal  and  spatial  filtering  affects  this 
distribution.  Maps  of  the  distribution  of  FSLE  for  a  fixed 
scale  are  often  used  to  identify  high-dispersion  regions  and 
mixing  boundaries.  In  particular,  the  spatial  distribution  of 
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Fig.  8.  Effect  of  space-smoothing  on  (a)  the  relative  dispersion  and  (b)  on 
the  FSLE  curve  for  six  different  values  of  a  The  trajectories  near  the 
coastline  have  been  discarded  to  avoid  spurious  dispersion,  and  the 
difference  in  D2  and  X  is  highlighted  by  the  new  (black)  and  original  (gray) 
full  field  curves. 

the  Lagrangian  structures  can  be  identified  by  plotting 
FSLE  maxima  from  forward  and  baekward-in-time  advee- 
tion  of  synthetic  drifters,  as  explained  in  d'Ovidio  et  al. 
(2004)  and  Moleard  et  al.  (2005).  To  this  end,  the  FSLE 
map  corresponding  to  the  model  simulation  on  September 
6,  2002  is  calculated  based  on  6  days  of  forward  and  back¬ 
ward  drifter  adveetion  using  an  initial  pair  distance  of 
0.45  km  and  y.  —  15.  This  selection  of  parameters  provides 
insight  on  the  structure  of  fine-scale  Lagrangian  features, 
while  allowing  them  time  to  fully  emerge. 

The  spatial  distribution  of  FSLE  for  the  full  model 
velocity  field  is  plotted  in  Fig.  9a.  The  map  clearly  shows 


the  inhomogeneity  in  the  dispersion  field  and  the  confine¬ 
ment  of  rapid  particle  pair  separations  to  a  few  distinct 
regions  and  isolated  flow  features.  The  vast  majority  of 
particle  pairs  undergoing  rapid  separation  to  scales  of 
order  15-30  km  initiate  in  confined  regions  along  the 
WAC,  in  the  Croatian  islands,  and  at  the  boundaries  of 
entraining  or  detraining  eddies  near  inter-gyre  boundaries 
and  behind  the  Gargano  Peninsula.  The  WAC  is  elearly  a 
controlling  influence  on  the  overall  relative  dispersion  in 
the  Adriatic.  The  intense  foliation  of  Lagrangian  structures 
shown  by  the  overlapping  of  FSLE  maxima  mark  the  com¬ 
plicated  pathways  by  whieh  Lagrangian  parcels  are 
entrained  and  detrained  from  the  boundary  eurrent.  In  that 
sense,  the  tangle  of  the  FSLE  curves  indicates  the  spatial 
extent  of  the  fluctuations  in  the  WAC  and  the  transport 
among  interacting,  small-scale  eddies  there.  Notably,  a 
vanishingly  small  number  of  model  particle  pairs  initiated 
in  the  gyre  interiors  separate  to  the  modest  a  —  15  seale 
in  the  6-day  integration  period.  The  number  of  identifiable 
‘eddy'  features  external  to  the  boundary  currents  as  shown 
by  isolated  transverse  intersections  of  forward  and  back¬ 
ward  in  time  FSLE  maxima  is  small.  Three  of  the  strongest 
are  located  in  the  southern  gyre  with  another  marked  by 
the  figure  eight  pattern  near  (200,625)  in  grid-point  units. 

The  main  effect  of  the  temporal  filtering  as  shown  in 
Fig.  9b-d,  is  to  remove  the  tangling  eaused  by  the  time- 
dependent  eddies,  filaments,  and  meandering  of  the 
WAC.  The  details  of  these  structures  are  considerably 
changed  for  smoothing  windows  as  short  as  1.5  days. 
The  other  isolated  eddy-mixing  regions  are  longer  lived 
and  persist  in  the  FSLE  maps  for  5-day  and  longer  averag¬ 
ing  times.  What  structure  persists  along  the  main  edges  of 
the  WAC  under  time  filtering  can  be  identified  with  distinet 
regions  of  high  strain  rate  in  the  mean  field,  for  example, 
the  mean  separation  bubble  behind  the  Gargano  peninsula 
and  the  flow  around  Isole  Tremiti  north  of  Gargano. 

Another  way  of  looking  at  this  is  by  plotting  the  prob¬ 
ability  density  distribution  (pdf)  of  i  at  different  separation 
distances.  Fig.  10  shows  that  at  distances  up  to  approxi¬ 
mately  20  km,  temporal  filtering  influences  the  pdf  of  t 
sueh  that  the  peak  taking  place  at  short  times  is  redueed 
and  the  pdf  moves  to  longer  time  scales  (Fig  10a).  In  con¬ 
trast,  temporal  filtering  has  no  impact  on  the  pdf  for  larger 
scales  (Fig.  10b).  Thus,  it  is  clear  that  while  the  complex 
fine-seale  instabilities,  eddies,  and  filaments  in  the  surface 
velocity  field  control  relative  dispersion  at  seales  smaller 
than  the  thickness  of  the  boundary  currents,  %20  km,  the 
shear  along  the  edges  of  the  boundary  currents  and  persis¬ 
tent  mean  flow  separation  features  entirely  dictate  the 
large-scale  dispersion.  As  sueh,  when  these  instabilities 
are  smoothed  out  by  time  averaging,  the  small-scale  disper¬ 
sion  is  affected  while  large-scale  relative  dispersion  remains 
intact. 

On  the  other  hand,  the  spatial  filtering  up  to  a  ^  5  km 
(Fig.  9e-g)  maintains  the  time-dependent  component  of 
the  tangling  of  the  Lagrangian  structures,  but  reduces  the 
intensity  of  local  velocity  gradients.  The  signatures  of  the 
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(b)  X(lp1  5d) 


(c)  X(lp5d) 


(d)  X(lp30d) 


Fig.  9.  Spatial  distribution  of  FSLEs  forward  6  days  in  time  starting  from  September  6,  2002  (positive)  and  backward  6  days  in  time  (negative)  calculated 
with  an  initial  pair  distance  of  0.45  km,  and  x  15  (units  are  h  ')  from  (a)  the  raw  field,  (b)  with  1.5-day  time  filtering,  (c)  with  5-day  time  filtering, 
(d)  with  30-days  time  filtering,  (e)  with  spatial  smoothing  using  a  2.5  km,  (0  a  5  km,  and  (g)  o  10  km 


WAC  and  corresponding  eddies  dissipate  for  a  =  10  km 
(Fig.  9f),  as  the  overall  sheer  is  reduced  and  spatially  dis¬ 
tinct  separation  regions  are  averaged  out  completely.  The 
pdf  of  i  (Fig.  10e)  behaves  in  a  similar  way  for  spatial  fil¬ 
tering  to  temporal  filtering  for  a  ^  10  km.  This  is  when 


small-scale  turbulent  structures  are  being  eliminated  from 
the  flow  field.  But  for  larger  a  values,  the  pdf  changes  sig¬ 
nificantly,  with  the  peak  moving  toward  longer  times  and 
decreasing  in  amplitude.  This  is  due  to  the  overall  reduc¬ 
tion  in  the  velocity  gradients  responsible  for  fast  disper- 
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days  days 


Fig.  10.  Comparison  of  the  distributions  of  dispersion  times  r,  defined  as  the  time  required  for  a  particle  pair  to  separate  from  8  to  for  the  full  and 
filtered  velocity  fields.  Two  separation  distances  are  chosen  to  highlight  the  filtering  effect  on  the  different  scales  of  the  flow:  (a)  2.1  2.5  km,  temporal 
filtering,  (b)  18.8  22.6  km,  temporal  filtering,  (c)  2  1  2  5  km,  spatial  filtering,  and  (d)  18.8  22.6  km,  spatial  filtering. 


sion.  The  data  indicate  that  this  occurs  for  averaging  scales 
smaller  than  the  width  of  the  boundary  currents  (Fig.  10d). 
Once  the  isolated  small-scale  spatial  structures  and  shear 
dispersion  zones  are  eliminated,  the  remaining  simpler  flow 
field  (lower  panel.  Fig.  7e)  combined  with  time  variability 
enforce  an  exponential  regime  in  relative  dispersion  that 
is  characteristic  of  chaotic  advection  (Fig.  8b). 

A  final  set  of  experiments  is  carried  out  to  further  differ¬ 
entiate  the  roles  played  by  persistent  mean  structures  and 
time-dependent  fluctuations.  This  is  done  by  examining 
the  relative  dispersion  statistics  resulting  solely  from  the 
time-dependent  component  of  the  flow  produced  by 
removing  the  monthly  mean  field.  This  fluctuating,  or  ‘tur¬ 
bulent’  field  is  given  by  u,(x,  /)  =  u(x,  /)  -  u(x),  where 
u(x,/)  is  the  monthly  mean  velocity. 

As  discussed,  low-pass  temporal  filtering  results  indicate 
that  the  large-scale  relative  dispersion  is  dominated  by 
mean  shear  zones  concentrated  in  the  boundary  currents. 
The  procedure  of  removing  the  mean  flow  is  therefore  sim¬ 
ilar  to  a  high-pass  filtering  in  this  case  (albeit  not  equiva¬ 
lent,  since  there  may  still  be  large  scale  residuals  from  the 
boundary  currents).  One  would  expect  the  results  to  be  sig¬ 
nificantly  different  with  respect  to  those  obtained  from  tem¬ 
poral  low-pass  filtering.  The  resulting  dispersion  and  FSLE 
curves  arc  plotted  in  Fig.  1 1,  and  display  behaviors  remark¬ 


ably  opposite  to  the  low-pass  experiments:  both  D2  and  X 
are  identical  to  the  full  flow  dispersion  at  small  temporal 
(/  <  2  days)  and  spatial  (<5  <  5  km)  scales,  while  the  disper¬ 
sion  at  larger  scales  changes  This  is  indeed  the  only  exper¬ 
iment  where  the  relative  dispersion  at  scales  smaller  than 
the  radius  of  deformation  remains  unchanged.  Therefore, 
the  results  indicate  that  the  shear  in  the  mean  flow  docs 
not  play  a  role  in  the  dispersion  at  small  scales.  After  20 
days,  D2  approaches  a  power  law  of  D2  ~  /*  \  and  a  diffu¬ 
sive  regime  after  1  month.  The  FSLE  at  scales  beyond 
20  km  displays  a  pure  ballistic  regime  X  l  0.  While  D2 
and  X  differ  regarding  the  dispersion  laws,  both  metrics 
show  that  when  persistent  mean  flow  structures  arc 
removed,  the  relative  dispersion  is  considerably  slowed  at 
large  temporal  and  spatial  scales.  The  distinct  ballistic 
regime  observed  in  the  FSLE  at  meso-  and  largcr-scalcs 
can  be  attributed  to  a  highly  organized  fluctuating  velocity 
field  characterized  by  long  correlation  times. 

5.  Summary  and  conclusions 

Quantification  of  relative  dispersion  is  important  for 
practical  problems  in  geophysical  fluid  dynamics  and  is  a 
classical  area  of  study  in  fluid  mechanics.  Several  regimes 
of  relative  dispersion  have  been  identified  on  a  theoretical 
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Fig.  11.  (a)  Relative  dispersion  D2U)  and  (b)  FSLE  A(r5)  curves  in  the  case 
of  detrended  velocity  field. 


basis  for  canonical  turbulent  Hows,  but  it  remains  unclear 
to  what  extent  these  regimes  arc  realized  in  the  ocean 
and,  in  particular,  in  complex  coastal  flows.  Investigation 
of  relative  dispersion  in  the  ocean  has  been  traditionally 
limited  by  the  limited  number  of  drifter  pairs  launched  with 
small  initial  separations.  Recent  advances  in  the  realism  of 
coastal  and  ocean  circulation  models  allow  investigation  of 
relative  dispersion  based  on  synthetic  drifters. 

Our  main  objective  is  to  study  relative  dispersion  in  the 
Adriatic  Sea  circulation  using  an  ocean  model  that 
employs  not  only  very  high  spatial  and  temporal  resolu¬ 
tion,  but  is  also  driven  by  realistic  surface,  tidal,  river, 
and  boundary  forcing.  The  Adriatic  is  a  scmi-cneloscd 


(marginal)  sea,  the  circulation  of  which  is  characterized 
by  strong  boundary  currents  and  interior  gyres.  The  sur¬ 
face  flows  in  the  Adriatic  Sea  arc  subject  to  strong  wind 
forcing  that  is  modulated  by  the  surrounding  land  topog¬ 
raphy.  The  buoyancy  input  by  the  Po  and  other  rivers 
leads  to  a  strong  WAC,  and  its  baroclinic  instabilities 
arc  known  to  facilitate  the  formation  of  eddies  and  fila¬ 
ments  that  should  presumably  be  very  important  for  rela¬ 
tive  dispersion  at  the  mcso-scalc  range.  As  such,  the 
Adriatic  is  an  interesting  and  challenging  place  to  explore 
this  problem.  Relative  dispersion  in  the  Adriatic  Sea  has 
been  subject  to  prior  investigation  by  Lacorata  ct  al. 
(2001)  using  37  real  drifters  and  an  idealized  kinematic 
numerical  model.  A  secondary  objective  of  this  study  is 
to  quantify  the  impact  of  the  uncertainties  that  can  arise 
from  numerous  errors  in  the  numerical  modeling  and/or 
forcing  data  on  the  resulting  Lagrangian  statistics.  This 
is  done  simply  by  applying  temporal  and  spatial  filtering 
to  the  model  Eulerian  velocity  fields. 

Relative  dispersion  based  on  the  raw  model  data  shows 
D^(r)~fy,  where  y  =  1.9  using  original  pairs  with 
Z)0  =  5.  lkm,  and  y  =  2.2  using  chance  pairs  with 
D0  =  1  km  for  t  ^  10  days  (a  different  type  of  launch  using 
chance  pairs  with  D0  =  1  km  yields  y  =  2.4),  while  the 
FSLE  shows  /  ~  S  0  8  (corresponding  to  D2(/)  ~  t2  5)  for 
S  ^  20  km.  As  such,  a  super-diffusive  relative  dispersion 
is  documented  at  these  scales  that  is  close  to  ballistic  based 
on  Dr(t)  and  between  the  ballistic  and  Richardson’s 
regimes  based  on  the  FSLE.  The  value  A  ~  t>~0,8  is  in  very 
good  agreement  with  results  from  Lacorata  ct  al.  (2001) 
based  on  real  drifters.  At  smaller  spatial  (and  faster  tempo¬ 
ral)  scales,  no  clear  dispersion  law  is  detected,  in  particular, 
an  exponential  regime  characteristic  of  chaotic  advection 
and/or  enstrophy  cascade  is  not  realized.  Analysis  of 
NCOM  output  from  additional  months  and  a  different  year 
of  simulation  show  no  significant  change  in  the  scaling 
laws,  which  indicates  that  the  flow  features  responsible 
for  the  overall  relative  dispersion  persist  across  seasons. 
Nevertheless,  it  should  be  noted  that  seasonality  is  com¬ 
puted  based  on  chance-pairs.  This  approach  undersamplcs 
the  regions  of  relatively  lower  activity  which  are  possibly 
more  sensitive  to  the  wind  regimes  and  thus  to  the  seasonal 
variability  of  the  forcings. 

Temporal  filtering  of  the  velocity  field  with  a  time  win¬ 
dow  up  to  30  days  only  affects  relative  dispersion  at  scales 
smaller  than  20  km  by  reducing  their  rate  of  dispersion  by 
half.  Relative  dispersion  at  scales  larger  than  20  km 
remains  intact.  In  contrast,  spatial  smoothing  with  a 
Gaussian  filter  with  a  up  to  20  km  affects  dispersion  at 
almost  all  scales  of  motion.  The  primary  difference  between 
temporal  and  spatial  smoothing  is  that  spatial  smoothing 
modifies  the  coherent  structures  of  the  mean  flow  field,  in 
particular,  the  shear  in  the  boundary  currents  for  smooth¬ 
ing  at  the  scale  of  these  currents.  As  such,  we  conclude  that 
large-scale  relative  dispersion  is  mainly  controlled  by  the 
mean  flow  in  the  Adriatic  Sea.  The  control  exerted  by  the 
large-scale  flow  on  the  tracer  evolution  is,  at  least,  not 
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inconsistent  with  the  ideas  put  forth  by  Bennett  (1984)  in 
that  the  slope  of  the  kinetic  energy  spectrum  is  steeper  than 
—  3  although  the  model  velocity  is  far  removed  from  any 
traditional  turbulence. 

Relative  dispersion  at  scales  smaller  than  20  km  and  10 
days  shows  no  clear  fit  to  the  traditional  ideas  and  both  the 
overall  statistics  and  spatial  structure  are  found  to  be  extre¬ 
mely  sensitive  to  the  resolution  and  sampling  of  the  WAC 
region.  When  spatial  smoothing  is  strong  enough  to 
weaken  the  influence  of  shear  dispersion  from  the  bound¬ 
ary  currents,  the  remaining  regime  appears  to  be  domi¬ 
nated  by  chaotic  advection.  This  implies  that  all  chaotic 
advection,  shear  dispersion,  and  small-scale  turbulent  fea¬ 
tures  such  as  filaments  and  eddies  play  a  role  in  relative  dis¬ 
persion  at  these  scales.  Since  these  scales  are  the  most 
important  for  environmental  problems  with  societal 
impact,  we  conclude  that  relative  dispersion  at  these  scales 
should  be  most  effectively  attacked  by  using  very  realistic, 
high-resolution  ocean  models.  The  accuracy  of  these  mod¬ 
els  needs  to  be  tested  and  improved  using  observations  at 
submcso-scalcs. 
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